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The standard model of cosmology ACDM assumes general relativity, flat space, and the presence 
of a positive cosmological constant. We relax these assumptions allowing spatial curvature, a time- 
dependent effective dark energy equation of state, as well as modifications of the Poisson equation 
for the lensing potential, and modifications of the growth of linear matter density perturbations 
in alternate combinations. Using six parameters characterizing these relations, we check ACDM 
for consistency utilizing cosmic microwave background anisotropics, cross correlations thereof with 
high-redshift galaxies through the integrated Sachs- Wolfe effect, the Hubble constant, supernovae, 
and baryon acoustic oscillation distances, as well as the relation between weak gravitational lensing 
and galaxy fiows. In all scenarios, we find consistency of the concordance model at the 95% confi- 
dence level. However, we emphasize that constraining supplementary background parameters and 
parametrizations of the growth of large-scale structure separately may lead to a priori exclusion of 
viable departures from the concordance model. 
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I. INTRODUCTION 

The detection of the late-time acceleration of our Uni- 
verse challenges the known laws of physics. General 
relativity and the standard model of particle physics can- 
not amount to the observed cosmic expansion unless we 
are willing to accept a seemingly random constant in the 
Einstein field equations or the presence of an unknown 
form of energy, along with large amounts of dark mat- 
ter. The concordance model assumes such a cosmologi- 
cal constant or vacuum energy and provides a simple but 
successful description for the observed Universe. Given 
the lack of a complete understanding, it is important to 
repeatedly test this model against observations with the 
ambition of distinguishing between different explanations 
of the observed cosmic acceleration, e.g., a cosmological 
constant, dynamical dark energy, or a modification of 
gravity, and rule out or constrain the different models. 
See [j-lal for recent reviews on tests of nonstandard cos- 
mologies. 

Phenomenological parametrizations for departures 
from the concordance model have been studied in [cl-fTsI . 
In this paper, we shall mainly adopt the parametriza- 
tion and notation of [l2|, [l^ . In addition to the usual 
cosmological parameters, including spatial curvature, we 
use five phenomenological parameters quantifying mod- 
ifications of the Poisson equation for the lensing poten- 
tial, modifications of the growth of linear matter density 
perturbations, as well as a time-dependent effective dark 
energy equation of state. The parameters are based on 
models of modified gravity, but may also describe prop- 
erties of dark energy and provide a framework to search 
for phenomena which may indicate new physical effects 
in current and future cosmological observations. 

We conduct a Markov chain Monte Carlo (MCMC) 
study of this parameter space using data from the cos- 
mic microwave background (CMB) anisotropics, super- 
novae distances, the baryon acoustic oscillations (BAO) 
distances, and the Hubble constant. We also utilize infor- 



mation from the cross correlation between high-redshift 
galaxies and the CMB through the integrated Sachs- 
Wolfe (ISW) effect, as well as a probe of the relation 
between weak gravitational lensing and galaxy flows. For 
the predictions of the CMB anisotropics, we connect our 
phenomenological parameters to theparametrized post- 
Friedmann (PPF) framework [13, Ey] and use its im- 
plementation into a standard Einstein-Boltzmann linear 
theory solver [T7| . 

In 211 we define the phenomenological parametriza- 
tion of the modifications to standard cosmology and ex- 
plain their implications on various cosmological probes 
in ^III Al Modifications to the ISWWLL O, [13 code 
used for the galaxy-ISW (gISW) cross correlation ob- 
servations are specified in i jlll A 41 For the connection 
of our parametrization to the PPF formalism, we refer 
to the appendix. We present the results of our MCMC 
study in illll Bl and discuss them in ijlVI 



II. PHENOMENOLOGICAL MODIFICATIONS 

We consider scalar linear perturbations of the 
Friedmann-Lemaitre-Robertson- Walker background in 
the longitudinal gauge, hence 



(1 + 2'^)dt^ + a^{l + 2^)dx^, 



(1) 



where diP is the unperturbed spatial line element with 
curvature k and c = 1. Motivated by modified grav- 
ity models, effects from a nonstandard cosmology may 
introduce the following three time- and scale-dependent 
phenomenological modifications on the background and 
at quasistatic scales of linear perturbations (see, e.g., [lol - 

EH nam; cf. i): 

• A deviation from the ACDM expansion history, 
parametrized by an effective dark energy equation 
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of state Wcs{a), 



H 



+(1 - a„ - f^k)a"^''"'="'^''^+^^ . (2) 

• An effective change of the Newton's constant, 
which we can parametrize by a free function 
Q{k,a), yielding a generahzed Poisson equation, 



(3) 



• A difference in the scalar linear potentials ^' and $ 
parametrized by the free function r](k, a), 



^' = -(1 + 77)$. 



(4) 



Deviations from Q = 1 and ij = 0, the standard values, do 
not necessary indicate modifications of general relativity 
but might, for instance, also originate from contributions 
of nonmatter fluids to the Poisson equation {Q ^ 1), e.g., 
from clustering of dark energy or interactions between 
the dark components (see, e.g., [l^), or nonvanishing 
anisotropic stress (77 ^ 0) (see, e.g., 0, H^). Dark en- 
ergy models other than the cosmological constant also 
predict departures from Wcs{a) = — 1. Note that a dark 
energy model can always be constructed to be formally 
equivalent to a modification of gravity through an effec- 
tive dark energy stress-energy tensor. The parameters 
introduced here cannot distinguish between the two de- 
scriptions. Consideration of microphysical aspects may, 
however, indicate which is the more reasonable picture 

{see, e.g., MM)- 
We can combine Eqs. ([3]) and ^ to obtain [1^ 



2a~ 



(5) 



where Y. = Q{1 + 77/2) and $_ = ($- *)/2. This re- 
lation is in particular sensitive to weak lensing measure- 
ments and the ISW effect can be used to probe its time 
evolution. Modifications as in Eqs. ^ and (|4]) lead to 
changes in the growth of the linear matter overdensity 
perturbation Am, which we parametrize via the growth 
index 7 



where ilm(i) 
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(6) 



We also introduce here 



the parameter ^ to account for growth rates beyond unity, 
given 7 > and < < 1, as can be observed in 
scalar-tensor theories~[l3l| . Hence, detection of ^ > may 
indicate the presence of an additional attractive force. 

We decide to use the quantities defined in Eqs. ([5]) 
and ([6]) to scan the cosmological data for departure from 
standard theory and replace the two free functions Q and 
77 by S and (7,0- Note that ^ should not be interpreted 
as a new parameter in the gravitational dynamics like Q 



and 77, but rather as the separation 7 — )■ 7 + P, where we 
assume the specific form P = ln(l -I- ^)/ ln57ni(a). 

From the ordinary differential equation, which de- 
scribes the correct behavior of the linear matter over- 
density perturbations in the quasistatic regime. 
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where T = Q(l + 77) = 2E(1 + 77)7(2 + 77) and primes 
denote derivatives with respect to In a, we derive 



^ = -r!m(a)27-i(i + ^)2 
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where we assumed constant 7 and ^. Given S and (7, 
we can use Eq. ([8]) to derive Q and rj and thus define 
our modifications to ACDM in the framework of Eqs. ([2]) 
through (gl). Note that 



(9) 



whenever a*^'""'^"^ » 1 and w^ff(hia) <C 1 as a ^ 1. 
If 1^1 ^ 1 at a ^ 1, standard gravity is reproduced. 
Otherwise, modifications to gravity persist up to high 
redshifts. 



A. Parametrization 

It is a difficult task to find general functions with a 
minimal set of free parameters that are flexible enough 
to capture the wealth of possible modifications in Eqs. ([2]) 
through (gl) or equivalently in Eqs. ©, ©, and ©. 
Rather than to construct such a function for each rela- 
tion, we decide to use a low number of five parameters, in 
addition to spatial curvature, and examine the combina- 
tional aspects of introducing modifications in each of the 
relations. We restrict to constant and time-dependent 
modifications and choose the parametrization in a way 
that ACDM is contained in the parameter space. 

For the expansion history we consider the parametriza- 
tion 



xa 



exp [3wa{a- 1)], 



(10) 



where the dark energy equation of state is given by 0, 01 



WDE(a) = u'o + (l-a)wa, (11) 
3 / ' l + ^DE(aO ^^,l 



a' 



(12) 



Note that this parametrization does not always provide 
a good fit to model predictions. We adopt the approach 
primarily due to its simplicity and wide usage. In the 
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limit (wQjWa) ^ (—1,0), the effective dark energy term 
in Eq. ([TU|) reduces to a cosmological constant. The 
parametrization is less successful, e.g., in the case of 
Dvali-Gabadadze-Porrati (DGP) gravity J25| . It provides 
a good approach for distance measures [Sl and the qua- 
sistatic regime for the self-accelerated branch, but it fails 
for the normal branch due to the appearance of a diver- 
gence in WDE{a)- Also note that choosing a specific form 
for wde may have a nontrivial effect on the constraints 
inferred for it (see, e.g., for a mo del- independent ap- 
proach) . 

Next, we need to define modifications to the rate of 
growth and the Poisson equation. We decide to use 
a constant growth index 70, which is a good approx- 
imation to general relativity, where 70 ~ 0.55, or the 



quasistatic regime of self-accelerating DGP gravity |27| . 
where 70 « 0.68 0, [2^, [2^. We further use a constant 
C = Coj which was found to provide a good fit to scalar- 
tensor theories where the scalar field couples to dark 
matter [isj . Equation ^ relates to the parametrization 
of [Tsj as d\ng^,/d\na = ^i7m(a)'*^, which was introduced 
to describe effects from early dark energy, modified grav- 
ity at high rcdshifts, or early acceleration. For these mod- 
els, g^, was found to be well-described by a constant p^ . 
Therefore, there is only limited correspondence to a con- 
stant ^ = Co = 0. For E, we use a parametrization that 
reduces to its general relativistic value S — 1 at early 
times. Hence, for our consistency test, we furthermore 
set 
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70, 
Co 

1 + Eoo- 



(13) 
(14) 
(15) 



For reference, we denote our nonstandard cosmologies by 
the extra parameters we allow to be free. Hence, when, 
e.g., taking fik and 79, as well as wq and Wq, to be free pa- 
rameters deviating from their standard values, we denote 
the according model by 7i/;k. Note that in the limit where 
{wo, wa, 70, Co, So) ^ {-1, 0, 0.55, 0, 0), Eqs. ©, 
([5]), and ^ reduce to general relativity with cold dark 
matter and a cosmological constant. From Eqs. ^ and 
PU)) . it is clear that when Wa ~ — wq, there is a de- 
generacy between Wa and Co, which manifests itself, in 
particular, at high redshifts. To explore the parameter 
space unclosed by free wq, Wa, and Co, we furthermore 
study three models where we fix 70 = 0.55 and Eq — 
while allowing the following degrees of freedom: 

(A) Co, W-DE = Wo + (1 - a)wa, 

(B) Co, W-DE = -1 + (1 - a)wa, 

(C) Co, u.DE = -l + Ao(l- 1.15a), 

such that ACDM is a limiting case in all of the three mod- 
els. The slope of ?i'DE(a) in model (C) is motivated by 
the best-fit values derived in mil Bl For numerical pre- 
dictions, we connect our parametrizations to the linear 
PPF framework as described in Appendix [Xj 



Note that we have ignored scale dependence of the pa- 
rameters, which is a good approximation within ACDM. 
Modifications of gravity such as f{R) gravity models [sol - 
[3^ . however, may introduce a strong scale dependence in 
these relations. The same holds for DGP gravity, where 
the deficiency is, however, restricted to near-horizon and 
superhorizon scales (3^. Although for these cases the 
parametrization is not very descriptive, it still serves as 
a useful tracer of nonstandard phenomenologies. Depar- 
ture from the standard parameter values may indicate in- 
consistencies in the concordance model and point toward 
new physical effects, which in turn have to be addressed 
with more developed theories. It has been pointed out, by 
conducting a principal component analysis (s^ , that such 
inconsistencies in Eqs. ^ through ^ arc more likely 
to be detected in scale-dependent modifications due to 
weaker sensitivity in the data to time-dependent devia- 
tions [ssl . [36j . A parametrization like the one described 
in this section is still capable of tracing trends in the ob- 
servations and disclosing new regions in the parameter 
space of possible modifications, especially when simulta- 
neously allowing for time-dependent modifications in the 
effective dark energy equation of state. 



III. CONSISTENCY CHECK 

We use a variety of cosmological data sets to check 
against nonstandard cosmology. First we use the CMB 
anisotropy data from the seven-year Wilkinson Mi- 
crowave Anisotropy Probe (WMAP) [33, the Arcminute 
Cosmology Bolometer Array Receiver (ACBAR) [38| . 
the Balloon Observations Of Millimetric Extragalac- 
tic Radiation and Geophysics (BOOMERanG) fiight in 
2003 (B03) [3§|, and the Cosmic Background Imager 
(CBI) [iOl. Next we employ data from the Supernova 
Cosmology Project (SCP) Union2 [4l| compilation, the 
measurement of the Hubble constant from the Super- 
novae and Ho for the Equation of State (SHOES) ^ 
program generalized by [4j] , and the BAO distance mea- 
surements of [31. Furthermore, we take gISW cross cor- 
relation observations using the iswwll code of [l^, [l^ , 
and the Eq measurement, probing the relation between 
weak gravitational lensing and galaxy flows, of [45| . 

Note that we restrict to data sets amenable to lin- 
ear perturbation theory. Nonlinear probes such as from 
the abundance of clusters and the full range of scales 
of weak gravitational lensing yield tight constraints on 
the parameters that quantify modifications to general 
relativity (see, e.g., j46l - l5l| ). In the case of /(i?) grav- 
ity models or DG P g ravity, nonlinear effects have been 
studied in [13, [53 - [60| and spherical collapse within var- 
ious dark energy models in, e.g., [6l|, [g^]- However, for 
phenomenological parametrizations like the one we use, 
nonlinear behavior has not been examined in full extent. 
Hence, applying standard scaling relations may lead to 
illusive conclusions. Furthermore, note that the inclu- 
sion of gISW cross correlations, along with the ISW ef- 
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feet in the CMB, and Eq, yield competitive results to 
constraints from nonlinear probes (cf., e.g., j63|). 

In mil Al we discuss the predictions for some of these 
observables for specific parameter values. In i ^IIIBI we 
present the results of a MCMC likelihood analysis, which 
is conducted with the publicly available COSMOMC [63 | 
package. 



A. Cosmological observables 

In this section we illustrate model predictions of vari- 
ous cosmological observables we use to derive our results. 
We chose the parameters of the various models that high- 
light results from the MCMC analysis. 

As our basic set, we choose a parametrization that sep- 
arates high-redshift and low-redshift constraints. Specif- 
ically we take six high-redshift parameters: the physical 
baryon and cold dark matter density fli,h^ and flch^, the 
ratio of sound horizon to angular diameter distance at 
recombination 0/100, the optical depth to reionization r, 
the scalar tilt n^, and amplitude As at = 0.002 Mpc~^. 
We extend this parameter set with alternate combina- 
tions of the six additional parameters, describing depar- 
tures from the concordance model: the spatial curvature 
density fJk, two parameters for the evolution of the effec- 
tive dark energy equation of state, wq and Wa, the growth 
index 70 and the scaling of the growth rate ^Oi cis well as 
the first-order of a time-dependent modification of the 
Poisson equation for the lensing potential Eq. 

We illustrate predictions from the maximum-likelihood 
ACDM model, as well as from the overall best-fit model 
(see ^IIIBp . Hereby, we derive the parameter values us- 
ing the full set of cosmological data. To demonstrate 
effects from the variation of a specific basic cosmological 
parameter on the observables, we use the corresponding 
ID-marginalized 68% confidence limits from the ACDM 
model while setting the complementary parameters to 
their best-fit values. For a supplementary parameter, we 
use its ID-marginalized 68% confidence boundary, when 
including it as the sole extra parameter while setting the 
basic parameters to their ACDM best-fit values. Note, 
however, that for the nonstandard cases, 70 is always si- 
multaneously varied with any of the extra parameters. 
Hereby, {wo,Wa) should be counted as one parameter; 
i.e., Wq and Wa are always either fixed to (—1, 0) or both 
considered free. 
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FIG. 1: Overall best-fit distance modulus with respect to 
the best-fit ACDM distance modulus. For illustration, the 
UNION2 data are binned into ten data bands with logarith- 
mic spacing. 
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FIG. 2: Deviation of the distance modulus with respect to 
the overall best-fit ACDM model. The dashed and dotted 
lines indicate upper and lower ID-marginalized 68% confi- 
dence limits of the according parameters, respectively, while 
other parameters are fixed to their ACDM best-fit values. 



tance defined by 



1. Distance measures 

The comparison of the magnitudes of high-redshift to 
low-redshift supernovae yields a relative distance mea- 
sure. Theoretical predictions for the distance modulus 
are related to the luminosity distance, dL{z) = (1 -I- 
z)r(z), where r(z) is the comoving angular diameter dis- 



sin[i/o\/=T\xW]/H"ovW, f^k<0, 

r(z) = <( x{z), a = 0, (16) 

sinh [HoV^xi^)] /Ho^A^l f^k > 0, 

where the comoving radial distance x is 

dz' 



Xiz) 



Hiz'Y 



(17) 
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FIG. 3: Inverse angular diameter distance dA (zed ~ 0.04) ~^ 
for the parameters qi £ {^lni,i^h,h,wo,Wa}- The solid seg- 
ments of the curves correspond to the ID-marginalized 68% 
confidence limits of the respective parameter. Only the ac- 
cording parameter is varied while other parameters are fixed 
to their best-fit ACDM values. The long-dashed line is the 
overall best-fit ACDM prediction and the dashed line corre- 
sponds to the prediction of the overall best-fit model. Note 
that due to the low effective redshift ZeB = 0.04, predic- 
tions for different values of Wa nearly overlap with the best-fit 
ACDM prediction. For fim and Jlk ID-marginalized 68% con- 
fidence intervals are too narrow to be distinguishable in the 
plot, i.e., fin, G (0.255,0.282) and G (-8.35 x 10"^2.10 x 
10~^), respectively. 



The supcrnovac magnitudes, once standardized, are re- 
lated to the distance by 

m = fi + M ^5 logio dL+M + 25, (18) 

where is in units of Mpc. The unknown absolute mag- 
nitude M of the supernovae is a nuisance parameter in 
the fit and is degenerate with Hq. Hence supernovae mea- 
sure relative distances within the set. In Fig. [1] we plot 
the prediction for the distance modulus from the overall 
best-fit model with respect to its counterpart from the 
best-fit ACDM model. The effect of varying parameters 
in the expansion history is illustrated in Fig. [2j 

The acoustic peaks in the CMB and the measure- 
ment of the local Hubble constant additionally provide 
absolute distance probes, which complement the rela- 
tive distance measure of the supernovae. For the Hub- 
ble constant, we utilize the SHOES measurement, Hq = 
74.2 ± 3.6 kms~^ Mpc""'^, which employs Cepheid mea- 
surements to link the low-redshift supernovae to the dis- 
tance scale established by the maser galaxy NGC 4258. 
In the analysis, we use the generalization of this measure- 
ment as a constraint on the inverse luminosity distance 
at zoff = 0.04, i.e., [H 

dAizesy^ = {I + Zcs) r {zcsy^ 

= {1 + Zcsf dL {Zcs)~^ 

~ (6.49 ±0.32) X 10-3 Mpc-\ (19) 

where dA is the angular diameter distance and the fidu- 
cial cosmology is h = 0.742, = 0.3, = 0, wq — —1, 
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FIG. 4: The ratio dz of the BAO measurement at redshifts 
2 = 0.2 (upper panel) and z = 0.35 (lower panel) for different 
parameter values. The horizontal dashed and long-dashed line 
corresponds to the overall best-fit model and best-fit ACDM 
prediction, respectively. They nearly overlap. The dot-dashed 
and dotted lines indicate upper and lower ID-marginalized 
68% confidence limits, respectively. These lines do not have 
any relevance in between the discrete parameter abscissae and 
only serve to indicate relative enhancement or suppression due 
to parameter variations with respect to the best-fit ACDM 
model. 



= 0. Figure 13] demonstrates predictions for (i^(2cff)~"'^ 
and effects on the observable from varying parameters in 
the expansion history. The gray band corresponds to the 
1(7 region of the SHOES measurement, where the solid 
gray line indicates the mean value. 

Further, we apply the BAO distance measurement 
of [3l that is obtained from analyzing the clusterin g of 
galaxies from the Sloan Digital Sky Survey (SDSS) [63 
and the 2-degree Field Galaxy Redshift Survey (2dF- 
GRS) [HBl, constraining the ratio 



dv{z) ~ (l + 2)2/3d^(^)2/32l/3 



(20) 



at z = 0.2 and z = 0.35. Here Vg {zd) denotes the co- 
moving sound horizon at the baryon drag epoch Zd- In 
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Fig.m we plot predictions for dz at z = 0.2 and z = 0.35 
and effects on the observables from varying parameters 
in the expansion history. The gray bands correspond to 
the la region of the BAO distance measurement, where 
the solid gray lines indicate the mean values. 

These probes place tight constraints on the background 
parameters in Eq. pU)) and help to identify sources for 
the growth of structure and break degeneracies. 



2. The cosmic microwave background 



SDSS [691, yielding Eq = 0.392±0.065 [451 at the redshift 
z = 0.32 by averaging over scales R = (10 — 50)h~^ Mpc. 

Figure [5] illustrates predictions for Eq for different 
parameter values of £ {f^m, fik, 7o, ^o, lUa, Sq, ^o}, 
when varying only one parameter at a time and fixing 
the others to their overall best-fit ACDM values. ID- 
marginalized 68% confidence limits are obtained from the 
MCMC analysis in ijlllBI The gray shaded band corre- 
sponds to the la region of the Eq measurement with the 
gray solid line being the mean value. 



The CMB probes the geometry of the background his- 
tory as well as the formation of large-scale structure. The 
latter manifests itself on the largest scales through the 
ISW effect from the evolution of the gravitational poten- 
tial. To predict these effects we connect our paramctriza- 
tion to the PFF formalism and utilize the PFF modifi- 
cations to CAMB [63] implemented in (T^]. We config- 
ure the FPF parameters as described in Appendix [A] 
This connection and the incorporation of the PFF for- 
malism into a standard Einstein-Boltzmann linear the- 
ory solver yields an efficient way to obtain predictions 
from our paramctrization for the CMB. It also prevents 
violations of energy-momentum conservation and avoids 
gauge artifacts in our results. 

In Fig. [SJ wc plot the CMB angular temperature 
anisotropy power spectrum for the overall best-fit model 
with respect to the prediction of the best-fit ACDM 
model. Effects from varying chain parameters are il- 
lustrated in Fig. [6l Figure [7] illustrates the CMB 
temperature-polarization (TE) cross power spectrum for 
ACDM and the overall best-fit model. The errors denote 
diagonal elements of the covariance matrix and include 
cosmic variance and instrumental noise. The overall best- 
fit model provides a better fit to all of the CMB data sets. 
This has to be attributed not only to the relative suppres- 
sion of the power spectra at large scales with respect to 
the best-fit ACDM model but also to deviations at large 
angular multipoles. 



3. Weak gravitational lensing and galaxy flows 

The relationship of weak gravitational lensing around 
galaxies to their large-scale velocities has been proposed 
as a smoking gun of gravity . The advantage of such 
a probe lies in its insensitivity to galaxy bias and initial 
matter fluctuations. The expectation value of the ratio 
of galaxy-galaxy to galaxy-velocity cross correlations of 
the same galaxies yields an estimator Eg 



111. We have 



Eg 



V2(<If - $) 

U a 111 a ^- 



a„(a)^ 1 + e' 



(21) 

Recently this quantity has been measured analyzing 
70 205 luminous red galaxies (LRGs) from the 



4- Galaxy-ISW cross correlations 

The correlation between galaxy number densities and 
the CMB anisotropics can be used to isolate the ISW ef- 
fect in the CMB and has proven to be a useful probe for 
constraining modifications to standard cosmology (see, 
e.g., [51I, m, [zO, ItiI) . We utilize the publicly avail- 
able ISWWLL code jl8l [l9j , where we turn off weak lens- 
ing contributions to the likelihood, focusing only on the 
gISW constraints. The galaxies in this probe are col- 
lected from the Two Micron All Sky Survey (2MASS) 
extended source catalog (XSC) [zllzl, the LRG samples 
and photometric quasars (QSO) of the SDSS (tJ, and the 
National Radio Astronomy Observatory (NRAO) Very 
Large Array (VLA) Sky Survey (NVSS) ^ and corre- 
lated with the five-year WMAP [zil CMB anisotropics. 
The resulting 42 data points of gISW cross correlations 
arc divided into nine galaxy sample bins j (2MASS0-3, 
LRGO-1, QSOO-1, and NVSS) based on flux (2MASS) or 
redshift (LRG and QSO). These data points are a selec- 
tion of multipole bins from all samples, where the selec- 
tion is based on the avoidance of nonlinearities and sys- 
tematic effects from dust extinction, galaxy foregrounds, 
the thermal Sunyaev-Zel'dovich effect, and point source 
contamination to affect the gISW cross correlations [l^ . 

We evaluate the gISW cross correlation in the Lim- 
ber and quasistatic approximation, as it is done in the 
ISWWLL code used for the data analysis. The gISW cross 
correlation in this approximation reads 



C 



(^+1/2)2 

D^G 

dz 



dzf,{z)H{z) 
+ 1/2^ 



X{z) 



(22) 



Here, P{k) is the matter power spectrum today. D is 
the linear density growth rate defined by Am(fc,z) = 
/\ra{k,Q)D / D\z=o and G = T.D[1 + z)/Y., is the linear 
potential growth rate, where = and Oj <C 1. 

In ACDM, the Limber approximation becomes accu- 
rate at the percent level for £ > 10 and drops approxi- 
mately as f' at higher £ (see, e.g., [77l - l79| ). This condi- 
tion is satisfied by about 90% of the total 42 data points 
that are used in the ISWWLL code. The approximation is 
also valid for, e.g., DGP and f{R) gravity l5l|,[7l|. Apart 
from the multipole £, the error depends also on the width 




FIG. 5: Ratio of the prediction of the CMB temperature anisotropy power spectrum of the overall best-fit model (dashed line) 
with respect to its best-fit ACDM counterpart. 
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FIG. 6: Ratio of the prediction of the ID-marginalized 68% confidence limits of cosmological parameters with respect to 
the overall best-fit ACDM prediction for the CMB temperature anistropies. Only the respective parameter is varied. Other 
parameters are fixed to their best-fit ACDM values. Upper limits are indicated by dashed lines, lower limits by dotted lines. 
Oscillatory behavior at the acoustic peaks indicate shifts of peak position. Note that the first acoustic peak is located at £ ~ 220 
and the second at / ~ 540. Various degeneracies between parameters and sets of parameters are clearly identifiable. 
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FIG. 7: CMB temperature-polarization cross power spectrum 
of the overall best-fit model (dashed line) and the best-fit 
ACDM counterpart (solid line). 
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FIG. 8: Eg a,t z = 0.32 for the parameters qi £ 
{Qm , Ok , 7o , Wo , Wa , So , ^0 } • The solid segments of the curves 
correspond to the ID-marginalized 68% confidence limits of 
the respective parameter. Only the according parameter is 
varied. Other parameters are fixed to their best-fit ACDM 
values. The long-dashed line is the overall best-fit ACDM pre- 
diction and the dashed line indicates the overall best fit. For 
Slk and ^0 ID-marginalized 68% confidence intervals are too 
narrow to be distinguishable in the plot, i.e., f2k G (—8.35 x 
10"^, 2.10 X 10"^) and G (-6.41 x 10"^ 8.66 x 10"^), re- 
spectively. 



of the redshift distribution, which changes only little with 
modifications. Given the large errors of the currently 
available data points at low £, we conclude that it is safe 
to apply the Limber approximation and furthermore that 
it is very useful since it is numerically faster than an exact 
integration. 

The function fj{z) relates the matter density to the 
observed projected galaxy overdensity with fj{z) = 
bj{z)Ilj{z) in the absence of magnification bias. 11^(2:) 
is the redshift distribution of the galaxies and the bias 
bj{z) is assumed independent of scale (cf. [t^]) but de- 
pendent on redshift. 



FIG. 9: gISW cross correlations for the NVSS sample. The 
solid line is the prediction from the best-fit ACDM model, 
whereas the dashed line corresponds to its counterpart from 
the overall best-fit model. 



The term in brackets in Eq. (j22p . parametrized by 
Eqs. ([S]) and ©, can be rewritten as 



D^G 

dz 



where 



a„(^r(i + e)]s + (i + z)^!>i?2 

dz 



D 



exp 



n^{zV{i + o 

z + l 



dz 



(23) 



(24) 



For SDSS quasars, the derivation of fj{z) involves the 
modified linear growth factor given through Eq. ([^ and 
magnification bias, for which we use a modification of the 
lensing window function of [Tsj 



W{z,z') 



3EiJ^(l 



2H{z) 



izf 



dlnr(z" 



dx{z") 



(25) 



Note that the gISW analysis uses photometric LRG 
samples, whereas the Eq measurement is based on spec- 
troscopic LRG samples that do not overlap in redshift. 
Furthermore, the gISW signal is dominated by large 
scales and most of the error is caused by sampling vari- 
ance and shot noise of galaxies. The error on Eq is dom- 
inated by uncertainties in lensing and redshift space dis- 
tortions and most of the signal comes from small scales 
around 10/i~^ Mpc. Therefore, we can safely neglect cor- 
relations between the gISW and Eq data sets. 

We plot predictions from the overall best-fit model and 
the best-fit ACDM model for the gISW cross correlation 
in the NVSS sample in Fig. |9] and illustrate effects of 
varying the chain parameters in Fig. 1101 



B. MCMC likelihood analysis 

We turn to the likelihood analysis of the cosmological 
parameter space. As our elementary set of parameters we 
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FIG. 10: Ratio of the prediction of the ID-marginalized 68% confidence limits of cosmological parameters with respect to the 
overall best-fit ACDM prediction for the gISW cross correlation in the NVSS sample. Other parameters are fixed to their 
best-fit ACDM values. Upper limits are indicated by dashed lines, lower limits by dotted lines. 



use C = {^i6/l^,^7c/^^,6',T, risjln [lO^'^A^] ,70}, where for 
the concordance model Cacdm ~ C n {70 ~ 0.55}. We 
implement the following flat priors: f2;,/i^ g (0.01,0.1), 
VLji^ e (0.045,0.99), Q € (0.5, lOVr e (0.01,0.8), n, e 
(0.5,1.5), ln[lOi%] e (2.7,4) jsl, and 70 G (-5,5). 
In addition we use all the sets of combinations between 
the phenomenological parameters and spatial curvature 
p^ e -P, where V = {IT, So, ^o, f^k} and W = {wo,Wa}, 
being free and taking on their ACDM values {wq = 
— l,Wa = 0, Eg = 0,^0 = 0,17k = 0}. We assign the fol- 
lowing flat priors to them: wq € (—5, 5), Wa G (—10, 10), 
So e (-5,5), Co G (-0.5,0.5), and G (-0.3,0.3). As 
starting centers for pi G V we use the ACDM values. For 
the (C) model, we use Ao G (—20, 20) with starting center 
Ao = 0.5. 

The COSMOMC package employs the Metropolis- 
Hastings algorithm [8l|, |82| for the sampling and the 
Gelman and Rubin statistic R [s^ for testing the con- 
vergence. We require R— 1 < 0.01 for our runs with two 
or fewer extra degrees of freedom and R — 1 < 0.02 for 
three or more extra degrees of freedom [95j . 

The chain samples are used to infer marginalized prob- 
abilities and mean likelihoods of the posterior. The 
marginalized distribution is obtained from projecting the 
samples to the reduced dimensions of a subspace, ignor- 
ing information about the goodness of fit and skewness of 
the distribution in the marginalized directions. Averag- 
ing the likelihood for each point of a subspace produces 
the mean likelihood. If the two curves do not overlap, 
the distribution is not Gaussian or the priors are not 
flat. It is important to note that in this case, marginal- 



ized probabilities may be amplified due to a larger pa- 
rameter volume rather than by a better fit to the data. 
Furthermore, for skew distributions, we quote our ID- 
marginalized constraints in terms of minimal credible in- 
tervals [s^l rather than by marginalized confidence lim- 
its. The former indicate where the tails occupy equal 
fractions of the probability distribution, whereas the lat- 
ter are constructed in such a way that any point inside 
of the interval has a larger posterior than a point outside 
of it, being the more meaningful selection in the case of 
skew or multimodal distributions. 

We summarize our results in Table U through HVl For 
comparison of the goodness of fit between the different 
models, we quote 

-2AlnA„ax = 2\n{C^^^'^/Cn,,^), (26) 
-2A{lnC)s = 2((ln£ACDM)s-(ln/:)s), (27) 
Ws = ('C)s/(/:acdm)s (28) 

in Table IIVI where £max is the maximum likelihood of 
a model and (•)s denotes the average over the Monte- 
Carlo samples. Note that MCMC sampling may not pro- 
vide very accurate best- fit estimates (see, e.g., |6j|). The 
best-fit values may also have a much higher likelihood 
than the mean, but simultaneously be confined to a very 
small region of the parameter space . We also give the 
mean likelihood of the samples {C)s, which corresponds 
to taking the posterior distribution as prior in calculating 
the evidence and in contrary to the maximum likelihood 
is a quantity that penalizes fine-tuning. If {C)s is greater 
than unity, this usually suggests that, on average, the ex- 
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Parameter 


Flat 


Nonflat 


70 
100 fik 


0.57 ±0.10 0.53 


0.59 ±0.10 0.60 
-0.32 ±0.53 -0.81 


70 

Wo 
Wa 

100 


0.59 ±0.11 0.56 
-1.10 ±0.16 -1.12 
(-0.19, 1.03) 0.44 


0.62 ±0.12 0.60 
-1.05 ±0.18 -1.18 
(-0.79,0.99) 0.58 
-0.54 ±0.77 -0.21 


70 
10 Eo 
100 Qi, 


0.57 ±0.12 0.62 
0.02 ±0.93 0.12 


0.58 ±0.11 0.60 
-0.03 ±0.95 0.24 
-0.30 ±0.58 -0.30 


70 
100 Co 
100 


0.58 ±0.10 0.57 
0.11 ±0.76 0.52 


0.60 ±0.10 0.56 
0.83 ± 1.08 0.47 
-0.73 ±0.75 -0.59 


70 

Wo 
Wa 

10 Eo 
100 Qi, 


0.58 ±0.12 0.55 
-1.11 ±0.15 -1.13 
(-0.99, 1.05) 0.48 

0.09 ±0.92 -0.08 


0.60 ±0.15 0.59 
-1.05 ±0.18 -1.18 
(-0.78,1.13) 0.48 
-0.14 ±0.94 -0.20 
-0.52 ±0.93 -0.25 


70 
10 Eo 
100 Co 
100 S^k 


0.57 ±0.12 0.51 
0.07 ±0.97 -0.27 
0.15 ±0.80 0.33 


0.60 ±0.11 0.56 
0.01 ±0.94 0.18 
0.79 ± 1.07 0.57 
-0.71 ±0.77 -0.47 


70 

Wo 
Wa 

100 ^0 
100 f2k 


0.60 ±0.11 0.67 
(-1.31,-1.03) -1.13 
(0.36,1.45) 1.17 
(-1.46,3.17) 10.06 


0.63 ±0.12 0.60 
(-1.30,-0.93) -1.29 
(-0.46,1.36) 1.29 
(-0.58,2.16) 3.72 
-0.64 ±0.93 0.23 


70 

Wo 
Wa 

10 Eo 
100 Co 
100 fik 


0.60 ±0.11 0.60 

(-1.29,-1.01) -1.15 
(0.37,1.43) 1.17 
0.27 ±0.99 -0.38 
(-1.29,3.14) 6.33 


0.61 ±0.13 0.66 
(-1.30,-0.94) -1.24 
(-0.37,1.42) 1.26 

0.06 ±0.96 0.50 
(-0.67,2.35) 7.19 
-0.54 ±1.02 -0.41 



TABLE I: Mean, standard deviations, and best-fit values 
for the extra degrees of freedom in the 7[i«] [E] [^] [k] models 
from using the WMAP7, ACBAR, B03, CBI, UNION2, BAO, 
SHOES, gISW, and Eg data sets. Values in brackets de- 
note the ID-marginalized 68% minimal credible intervals [s^ 
quoted for parameter directions in the posterior distribution 
with distinctive skewness. 



tra parameters improve the goodness of fit to the data. 
This, however, has to be interpreted only rule of 
thumb (see, e.g., [i^). We only give constraints on the 
parameters Ci G C\{jo} for the 7W^, (A), (B), and (C) 
models in addition to ACDM (see Tables HI through Hill) . 
These models, except for the (B) model, exhibit an im- 
proved goodness of fit in terms of maximum and averaged 
likelihoods over the ACDM model. We numerically eval- 
uate Eqs. ([26| through ([28]) and quote the numbers for 
the -/w[Y.]C[k], (A), and (C) models in Table IVl For all 
other models, {C)s is smaller than unity, — 2A(ln£)s is 
positive, and — 2Aln£inax > —0.5. 



Under the assumptions made in ^11 Al our results show 
that constraints on the growth index parameter 70 are 
only weakly affected by the introduction of the other ex- 
tra degrees of freedom with consistency of ACDM, i.e., 
7o ~ 0.55, at the 68% confidence level. The constraints 
on 7o are dominated by the gISW data and are com- 
petitive to existing results in the literature derived from 
linear and nonlinear probes (cf., e.g., [46, 48, 63, 84]). 

In every scenario studied here, the standard values of 
all extra parameters, corresponding to the concordance 
model, lie within their ID-marginalizcd 95% confidence 
limits. Thus, we conclude that the ACDM model is con- 
sistent with the joint set of WMAP7, ACBAR, BOS, CBI, 
UNI0N2, BAO, SHOES, gISW, and Eq data under the 
assumption of the existence of extra degrees of freedom 
of the kind described in mi Al and form of priors given 
above. Note that constraints on a parameter direction 
depend a lot on the prior assumed for it or for the pa- 
rameters used to derive it (sec, e.g., ^^). Identifying 
adequate priors on extra degrees of freedom in the gravi- 
tational dynamics based on theoretical contemplations is 
the object of current research (see, e.g., [2^ Issj). 

For the maximum-likelihood ACDM model we obtain 
a fit of -21n£^,^DM ^ §099.3 using all of the data de- 
scribed above. The overall best-fit to the data was ob- 
tained in the chains of model (C) with — 2Aln£max = 
-7.3 with respect to the best-fit ACDM model. Note 
that in principle the best-fit values of model (C) are also 
contained in the parameter spaces of the 7w[S]^[k] and 
(A) models. The fact that none of these models accessed 
the best-fit point can be attributed to sampling errors in 
the chains. However, for the maximum-likelihood values 
of the 7w[I]]^[k] and (A) models, we also obtain increases 
in fit of - 2 A In /Imax < -2 over the best-fit ACDM model 
(see Table lIVp . In terms of p values this translates to 
3%, 18%, and 27% for consistency with ACDM in the 
(C), (A), and jwS^ scenarios, respectively, where for all 
other models this number is greater than 50%. How- 
ever, several best-fit parameter values of the 7U'[S]^[k], 
(A), and (C) models lie beyond their corresponding ID- 
marginahzed 68% confidence limits. In the case of, e.g., 
the jw^ model with parameter space CU{wo, Wa,^o}j the 
best fit of ^0 lies even beyond the ID-marginalized 99% 
confidence level, indicating that the best-fit model occu- 
pies only a very small subspacc of the parameter space. 

Remarkably, this increase in maximum likelihood 
(Axmax ^ ~2) seems to appear only in cases where we 
allow a free, positive, ^Qj phantom crossing in w-de, and 
wq « —Wa, i.e., an effective dark energy equation of 
state that drives toward a matterlike equation of state 
at early times. We further observe that in these sce- 
narios, there is no signature of 7^ 1. Quite to the 
contrary, the overall best fit obtained in model (C), i.e., 
Us = 0.996, portrays an almost perfect scale-invariant 
Harrison-Zel'dovich spectrum [Hjls^l- Moreover, we ob- 
serve a slight preference for smaller dark matter densities, 
smaller amplitudes for the primordial superhorizon power 
of the curvature perturbation, and larger ratios of sound 
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Parameter 


ACDM 




(A) 


100 Qth'^ 


2.229 ±0.050 


2.244 


2.237 ±0.056 


2.246 


2.231 ±0.056 


2.196 




0.1114 ±0.0030 


0.1121 


0.1066 ±0.0075 


0.0995 


0.1058 ±0.0077 


0.0935 


6 


1.0397 ±0.0024 


1.0406 


1.0443 ± 0.0063 


1.0510 


1.0453 ± 0.0068 


1.0560 


T 


0.085 ±0.014 


0.087 


0.086 ±0.015 


0.090 


0.086 ±0.015 


0.091 


Us 


0.960 ±0.012 


0.962 


(0.947,0.991) 


1.016 


(0.948, 0.992) 


1.011 


In [10^«A.] 


3.199 ±0.036 


3.199 


(3.034, 3.270) 


2.810 


(3.009,3.271) 


2.744 


70 






0.60 ±0.11 


0.67 






wo 






(-1.31,-1.03) 


-1.13 


(-1.29,-1.02) 


-1.07 








(0.36, 1.45) 


1.17 


(0.38, 1.47) 


1.09 


100 Co 






(-1.46,3.17) 


10.06 


(-1.42,3.63) 


9.67 




0.268 ±0.014 


0.269 


0.258 ±0.018 


0.233 


0.256 ±0.019 


0.224 


0"8 


0.806 ± 0.021 


0.812 


0.779 ±0.077 


0.749 


0.800 ± 0.063 


0.769 


Ho 


70.6 ± 1.3 


70.7 


70.7 ± 1.7 


72.3 


70.7 ± 1.8 


71.8 



TABLE II: The same as in Tabled but for the ACDM, "fw^, and (A) models and including the basic cosmological parameters 
along with matter density, power spectrum normalization aa, and the Hubble constant. 



Parameter 


(B) 


(C) 


(C)* 


100 Qbh'^ 

6 
r 

Us 

In [10^'^ As] 

Wa 

100 Co 
Ac 


2.236 ± 0.056 2.230 
0.1097 ±0.0059 0.1150 
1.0413 ±0.0048 1.0395 
0.085 ±0.014 0.086 
0.963 ±0.014 0.959 
3.179 ±0.059 3.221 
(-0.23,0.53) -0.20 
(-0.81,1.41) -0.18 


2.229 ±0.056 2.189 
0.1041 ±0.0073 0.0976 
(1.0474,1.0559) 1.0550 
0.086 ±0.015 0.081 
(0.953, 1.002) 0.996 
(2.830,2.958) 2.841 

(5.07, 7.64) 7.88 
(0.31,1.22) 1.02 


2.232 ±0.055 2.201 
0.1023 ±0.0067 0.0986 
1.0489 ± 0.0066 1.0544 
0.087 ±0.014 0.085 
(0.970, 1.018) 0.998 
(2.775, 2.982) 2.864 

(4.19,8.59) 7.68 
(0.71,1.25) 1.02 


era 
Ho 


0.267 ±0.015 0.274 
0.809 ± 0.060 0.829 
70.4 ± 1.7 70.9 


0.253 ±0.021 0.227 
0.811 ±0.062 0.818 
70.7 ± 1.8 72.5 


0.248 ±0.021 0.232 
0.788 ± 0.062 0.805 
71.0 ±2.0 72.0 



TABLE III: The same as in Table IH but for the (B) and (C) models, where (C)* denotes the case where no gISW data are 
used. The best-fit values quoted here for the (C) model define the overall highest maximum-likelihood model found for the 
combination of all of the data sets studied in this paper and are used to derive predictions for illustrations in i )III Al For the 
(C) models, the ID-marginalized minimal credible intervals correspond to the mode with larger mean likelihood. 



horizon to angular diameter distance at recombination 
with respect to parameter values inferred for the concor- 
dance model. If we additionally assume a flat Universe 
and matter density perturbations that relate to lensing 
potentials as i n AC DM, i.e., Sq = 0, these scenarios in- 
duce values for (£)s beyond unity, suggesting an increase 
of the goodness of fit through the inclusion of their extra 
parameters. 

The increase in maximum likelihood observed for mod- 
els with ^0 > and w — wq > 1 is mainly attributed 
to slightly better fits of the CMB anisotropy data on all 
scales and to a smaller extent also to marginally better 
fits in distance measures, as can be perceived from the 
figures in ijlll Al Merely the gISW data counteracts this 
trend in favor of ACDM and when removed from the 
joint data set for, e.g., the (C) model, here denoted by 
(C)* , we obtain the most extreme values for our statistics. 



Eqs. ^ through ^ (see Table |TVl), while ACDM pa- 
rameter values remain consistent at the 95% confidence 
level (cf. Table mi). 

FiguresfTTjandfT^show the marginalized likelihoods for 
the extra parameters and the parameters that exhibit a 
distinctive skewness or multiple modes in the case of the 
and (C) model, respectively. In Fig. [T21 we plot 2D- 
marginalized contours of the extra parameters wq, Wa, 
7o, and within the 710^ model. We also indicate the 
ID-marginalized minimal credible intervals on Aq and ^0 
obtained in model (C). Note that the best- fit points lie 
outside the 2D-marginalized 68% confidence level con- 
tours of ^0 and 7o, as well as of ^0 and Wa, which is due 
to the fact that the best-fit model occupies only a small 
parameter subspace and a distinctive skewness of the pos- 
terior distribution in the parameter directions. Note that 
there is a strong degeneracy between the parameters ^0 
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Statistic 


ACDM jw^ jw£,k 7wE^ 'yw'E^k (A) (C) Others 


(C)* 


— 2Aln £max 

-2A(In C), 


-5.2 -2.2 -3.8 -3.4 -4.9 -7.3 > -0.5 

0.9 2.9 2.2 3.4 0.0 -0.9 > 

1 2.6 0.6 1.0 0.6 2.3 6.5 < 1 


-6.8 
-2.9 
11.3 



TABLE IV: Comparison of the goodness of fit of the models with respect to ACDM, where InC^S^^^ = -4049.6, (ln£ACDM>s = 
—4052.8, and ln(£ACDM)s ~ —4051.9. (£)s > 1 suggests that the adoption of the extra parameters might be justified. (C)* 
denotes the case where no gISW data are used. 



and Wa as we expected from Eq. Figure [H] illustrates 
the overall best-fit parameter values from model (C) and 
jw^ for the modifications in the dark energy equation of 
state WDE(a), the factor J- that enters the ordinary dif- 
ferential equation of quasistatic matter overdensity per- 
turbations, the scaling of the effective Newton's constant 
Q, and rj that modifies the ratio of the scalar potentials. 
Also shown are the standard values for these quantities. 

We find consistency of the ACDM values at the 95% 
confidence level for all extra parameters in all scenarios. 
We, therefore, do not claim detection of nonstandard cos- 
mological effects or the necessity of extra degrees of free- 
dom based on the increases in the maximum and mean 
likelihoods (see Table ITV)) . This would require a physi- 
cally better-motivated model with well-understood priors 
and a more sophisticated analysis of the goodness of fit 
with, e.g., efficient exploration of the different modes and 
determination of the Bayes factor for model comparison 
(see, e.g., (sst - fQlj ). Note that even for a hypothetical, 
theoretically well-motivated model that would naturally 
access the best-fit region in the parameter space, there 
would only be a moderate preference over the concor- 
dance model, i.e., Ax^ax — It is not clear if such 
a model exists or may be developed, but we emphasize 
that we may a priori be excluding viable departures from 
the concordance model when constraining supplementary 
parameters for the background and growth of structure 
separately. 



95% confidence level. For specific combinations of the 
supplementary free parameters, we, however, find an in- 
crease in the maximum likelihood of up to Ax^^^ ~ — 7 
over the maximum-likelihood ACDM model in the joint 
data. We observe that this increase in maximum likeli- 
hood only appears in cases where we allow a free, posi- 
tive, ^0, phantom crossing in wdeCo), and wq w —Wa, i.e., 
an effective dark energy equation of state that drives to- 
ward a matterlike equation of state at early times. More- 
over, we find that in these scenarios there is no preference 
for rig ^ 1. For scenarios where we adopt the assumptions 
of a flat Universe and a standard relation of matter den- 
sity perturbations to lensing potentials, both maximum 
and mean likelihoods are greater than in ACDM, with 
a ma ximal ratio of mean likelihoods of the samples of 
(C)s — 6, suggesting that the introduction of their extra 
parameters might be justified. The increase in maximum 
likelihood can be attributed to better fits of the CMB 
anisotropy data on all scales and to a smaller extent also 
to better fits in distance measures. We therefore expect 
future CMB data from, e.g., the Planck mission [92| to 
yield more decisive constraints on our modifications (cf., 
e.g., ^^)- The gISW data sets counteract this trend in 
favor of the ACD M mo del and when removed, we obtain 
AXmax - -7 and {C)s ^ 11 for the (C) model. 



IV. DISCUSSION 

We allowed departures from the concordance model by 
introducing six additional degrees of freedom parametriz- 
ing a time-dependent effective dark energy equation of 
state, modifications of the Poisson equation for the lens- 
ing potential, and modifications of the growth of linear 
matter density perturbations, as well as spatial curva- 
ture. We constrained alternate combinations of these 
parameters by performing MCMC likelihood analyses on 
cosmological data amenable to linear perturbation the- 
ory. In particular, we utilized all of the CMB data, in- 
cluding the lowest multipolcs, its correlation with galax- 
ies and the comparison of weak lensing to large-scale 
velocities in addition to geometrical probes from super- 
novae and BAO distances, as well as from the Hubble 
constant. 

We find consistency of the concordance model at the 



Given the consistency of ACDM parameter values with 
their marginalized constraints, also under removal of 
gISW data, and the lack of a better-motivated theory 
and priors, we do, however, not claim detection of non- 
standard cosmological effects or the necessity of the in- 
troduction of extra degrees of freedom, but we emphasize 
that when constraining the background parameters and 
growth parameters separately, we may a priori be exclud- 
ing viable departures from ACDM. We leave the search 
for a well-motivated model that matches the above re- 
quirements and the analysis of its performance in con- 
frontation with new data, as well as the study of effects 
from scale-dependent deviations from the concordance 
model on our results to future work. We also point out 
that more sophisticated sampling methods may signif- 
icantly improve our statistical analysis and offer more 
information through efficient exploration of the differ- 
ent modes and the determination of the Bayes factor for 
model comparison. 
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wq vy„ fo ln[10"'A,] n, 

FIG. 11: Marginalized likelihoods for different parameters of the 710^ model that exhibit skew distributions. 




Ao fo ln[10'"A,] «, 

FIG. 12: The same as Fig. Illl but for the (C) model. The posterior distribution is bimodal. 
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Appendix A: Connecting our parametrization to the 
linear PPF framework 

Given the expansion history, the PPF framework [T^ . 
fl6j is defined by three functions and one parameter. 
From these quantities, the dynamics are determined by 
conservation of energy and momentum and the Bianchi 
identities. The defining quantities are g{a,k), which 
quantifies the effective anisotropic stress of the modifica- 
tions and distinguishes the two gravitational potentials, 
f({a), which describes the relationship between the mat- 



ter and the metric on superhorizon scales, and /gIi), 
which defines it in the linearized Newtonian regime. The 
additional parameter cr is the transition scale between 
the superhorizon and Newtonian behaviors. 

From the relations defined in i jTTl we infer g = —Tj/{rj + 
2) = 1 — and /g = (1 — in the quasistatic 

regime. At superhorizon scales, we impose f( = and 
we set the transition at the horizon, cr = 1. Given that 
g — >■ sufficiently fast for a < 1, if |^| ^ 1 and E — > 1, 
we apply the relation g = 1 — T/H to all scales. 

We implement these relations and the definitions in f|TT] 
into the PPF modified CAMB code This proce- 

dure produces the correct power spectra for ACDM and 
self-accelerating DGP gravity. The latter, however, only 
when restricting to its subhorizon effects. In connecting 
our parametrization to the PPF linear theory, we can 
take advantage of a fully consistent framework for mod- 
ifications in the gravitational dynamics. In this way, we 
can also prevent implicit violations of the conservation of 
energy and momentum and avoid gauge artifacts in our 
parametrization. 
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